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The  3D  compact  model  of  a  thermoelectric  cooler  in  the  CFD  approach  is  an  essential  technique  for  com¬ 
putationally  realistic  design  of  many  industrial  thermal  management  applications.  The  first  objective  of 
this  paper  is  to  develop  a  “black  box”-like  model  for  thermoelectric  coolers  as  well  as  a  method  to  for¬ 
mulate  the  effective  material  properties  for  the  compact  model  based  on  the  manufactures’  datasheets. 
Second,  a  highly  detailed  and  physical  thermoelectric  cooler  model  is  implemented  to  carry  out  virtual 
experiments  when  the  cooler  datasheet  is  not  available.  A  number  of  close  comparisons  validate  that 
the  compact  model  and  parameter  extraction  are  accurate  enough,  and  the  computation  time  and  size 
are  significantly  lower  than  the  physical  model.  The  compact  model  is  also  demonstrated  to  be  helpful 
in  evaluating  the  thermal  system  components  such  as  an  air-cooled  heat  sink,  so  that  the  system-level 
thermal  management  optimization  is  possible  without  complex,  multi-scale  computation. 
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1.  Introduction 

Thermoelectric  cooling  technology  can  reduce  the  temperatures 
of  electronic  module  packages  below  the  ambient  environment 
temperature  without  the  use  of  pumped  liquid  cooling  or  vapor- 
compression  refrigeration.  The  solid  state  refrigerators  are  hence 
particularly  appropriate  to  the  growing  situation  of  cooling  where 
reliability,  volume  and  complexity  are  of  prior  consideration  to  the 
cooling  coefficient  of  performance  (COP).  In  addition,  the  stringent 
temperature  control  required  by  some  optoelectronic  components 
can  also  be  provided  by  thermoelectric  coolers  (TECs),  the  cooling 
of  which  is  readily  controlled  through  adjusting  the  current  ap¬ 
plied.  A  constant  temperature  or  heat  flux  of  the  TEC  can  be  main¬ 
tained  by  the  tight  thermal  management  even  if  the  ambient 
thermal  condition  around  it  varies  dynamically.  For  these  benefits, 
there  has  been  an  increased  interest  in  the  application  of  TECs  to 
electronic  and  photonic  device  cooling  to  handle  the  heat  from 
the  high  power  density  devices,  and  TECs  have  actually  been  used 
in  many  diverse  applications  for  extending  the  ability  of  air/liquid- 
cooled  heat  sinks  [1-8]. 

The  successful  design  of  a  thermoelectric  cooling  system  neces¬ 
sitates  an  excellent  TEC  simulation,  which  ideally  should  be  able  to 
analyze  the  cooling  performance  and  predict  the  thermal  behavior 
with  sufficient  accuracy.  A  number  of  three-dimensional  (3D)  TEC 
models  have  been  reported,  in  which  the  geometric  detailed  repre¬ 
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sentation  of  the  entire  array  of  thermoelements  is  exactly  modeled 
[9-13].  However,  these  TEC  models  are  not  developed  within  com¬ 
putational  fluid  dynamics  (CFD)  software  packages,  and  hence  dif¬ 
ficult  to  be  integrated  into  a  larger  system-level  analysis  against 
the  fluidic  heat  exchange.  Also  based  on  the  true  solid  representa¬ 
tion  of  the  thermoelectric  device,  other  3D  models  have  been 
implemented  by  the  original  approach  for  a  CFD  simulation,  mak¬ 
ing  the  co-simulation  of  thermoelectric  devices  and  fluid  heat 
sources  convenient  [14,15].  In  [15],  the  popular  CFD  code  FLUENT 
is  numerically  equipped  with  a  thermoelectric  model  for  genera¬ 
tion  applications,  where  the  temperature  dependence  of  all  proper¬ 
ties  of  the  p-  and  n-type  materials,  the  thermo-electric  coupled 
multiphysics  field,  and  the  complete  thermoelectric  processes  of 
Seebeck,  Peltier,  and  Thomson  effects  integrated  with  Joule  source 
terms,  are  all  taken  into  account. 

The  integration  of  thermoelectric  devices  into  CFD  simulation 
environments  is  also  of  interest.  Here  a  major  interest  and  concern 
from  industry  is  the  practicality  and  cost  of  modeling.  While  the 
multiphysics  packages  are  ideal  for  modeling  each  n-type  and  p- 
type  element,  thermal  and  electrical  contact  individually  [9-15], 
the  cell/grid  size  should  be  on  the  order  of  microns  to  accurate  cal¬ 
culate  the  thermoelectric  device.  The  ideal  cell/grid  size  for  the  sys¬ 
tem  is  typically  orders  of  magnitude  larger  introducing  “grid 
mismatch”  making  integration  of  the  multiphysics  model  of  a  TE 
with  a  system  model  very  computationally  expensive  and  compli¬ 
cated  to  perform. 

Modeling  the  TE  module  as  a  single  bulk  block  or  “black  box”  to 
represent  the  TEC  instead  of  modeling  each  thermoelement 
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Nomenclature 

A 

uniform  cross-sectional  area  of  the  entire  TEC,  m2 

i 

n-type  material  figure  of  merit 

Ap 

An 

uniform  cross-sectional  area  of  p- type  leg,  m2 
uniform  cross-sectional  area  of  n-type  leg,  m2 

Z 

effective  device  figure  of  merit 

CFD 

computational  fluid  dynamics 

Greek  symbols 

COP 

coefficient  of  performance 

a 

effective  Seebeck  coefficient  of  the  compact  TEC  model, 

/ 

packing  fraction  of  total  TEC  area  covered  by  thermo- 

VIC1 

electric  elements 

ocp 

Seebeck  coefficient  of  the  p- type  materials,  V  K1 

i 

reduced  current 

GCn 

Seebeck  coefficient  of  the  n-type  materials,  V  K-1 

I 

electrical  current,  A 

AT 

=Th  -  Tc,  external  temperature  difference  between  the 

fm  ax 

the  electrical  current  in  the  device  that  gives  ATmax,  A 

hot  and  cold  sides,  K 

I< 

effective  thermal  conductance  of  the  TEC,  W  K_1 

A  Tmax 

the  maximum  temperature  difference  when  there  is  no 

l<, 

effective  parallel  thermal  conductance  within  the  TEC  of 

heat  absorbed  at  the  cold  side,  K 

one  pair  of  thermoelements,  W  K-1 

ATte 

=Thi  -  Td,  internal  temperature  difference  between  the 

l 

thermoelectric  element  length  or  height,  m 

hot  and  cold  junctions,  K 

L 

height  of  the  macro  block  (/  plus  the  height  of  the  top 
and  bottom  metal  interconnects),  m 

K 

effective  thermal  conductivity  of  the  compact  TEC  mod¬ 
el,  Wm^r1 

N 

number  of  couples  in  the  TEC  module 

Kt 

effective  parallel  thermal  conductivity  within  the  TEC  of 

Qjoule 

the  Joule  heating  density,  W/m3 

one  pair  of  thermoelements,  W  m_1 1<-1 

Qpelt 

the  Peltier  heat  generation  rate,  W/m3 

Kp 

thermal  conductivity  of  the  p- type  semiconductor 

Q 

the  rate  of  heat  transfer  at  the  side  being  cooled,  W 

material,  W  m-1  K-1 

Qmax 

the  maximum  rate  of  heat  transfer  at  the  side  being 
cooled  with  no  temperature  difference,  W 

Kn 

thermal  conductivity  of  the  n-type  semiconductor 
material,  W  m-1  K-1 

Qc 

total  rate  of  heat  transfer  from  the  heat  source,  W 

P 

effective  electrical  resistivity  of  the  compact  TEC  model, 

Qh 

total  rate  of  heat  transfer  to  the  heat  sink,  W 

Q  m 

R 

electrical  resistance  of  the  TEC,  C2 

Pi 

effective  parasite  electrical  resistivity  of  one  pair  of 

Ri 

effective  parasite  resistance  of  one  pair  of  thermoele¬ 

thermoelements,  Q  m 

ments,  C2 

Pp 

electrical  resistivity  of  the  p- type  semiconductor  mate¬ 

s 

=2na,  effective  lumped  Seebeck  coefficient  of  the  TEC 

rial,  Q  m 

black  box,  V  I<-1 

Pn 

electrical  resistivity  of  the  n-type  semiconductor  mate¬ 

T 

temperature,  K 

rial,  Q  m 

Tc 

temperature  of  the  TEC  cold  side,  K 

yP 

=  ^r,  relative  area  for  the  p- type  thermoelectric  ele¬ 

Th 

temperature  of  the  TEC  hot  side,  K 

ments 

Td 

cold  junction  temperature,  K 

fn 

=  relative  area  for  the  n-type  thermoelectric  ele- 

Thi 

hot  junction  temperature,  K 

mentsP 

TEC 

thermoelectric  cooler 

0C 

cold  substrate,  W_1  K 

V 

voltage  of  external  power  supply,  V 
p- type  material  figure  of  merit 

0n 

hot  substrate,  W-1  K 

individually  can  make  the  computation  much  easier  [16-23].  This 
kind  of  compact  models  is  able  to  cope  with  the  multiscale  issue  at 
the  interface  between  the  very  fine  mesh  of  the  thermoelement  re¬ 
gion  and  relatively  coarse  mesh  of  the  package  and  fluid  region.  Be¬ 
cause  the  TEC  module  as  a  whole  rather  than  each  pellet  is  modeled, 
the  TEC  gridding  can  be  significantly  simplified  and  the  grid  mis¬ 
match  can  thus  be  mitigated.  From  a  computational  efficiency 
standpoint,  full-field  CFD  and  heat  transfer  simulations  are  now  less 
expensive  but  more  feasible  as  the  number  of  volumes  as  well  as  the 
algorithm  complexity  is  reduced.  In  applying  compact  TEC  models, 
however,  retaining  the  accuracy  in  capturing  the  cooling  perfor¬ 
mance  is  as  important  as  the  simulation  simplicity. 

It  is  noteworthy  that  the  model  parameters  have  decisive  influ¬ 
ence  on  the  simulation  results.  Thus,  the  process  of  parameter  pro¬ 
curement  is  inextricable  to  the  method  of  building  a  compact 
model  itself,  and  actually  a  key  design  step  in  the  model  construc¬ 
tion.  The  critical  technique  here  is  how  to  determine  the  effective 
values  of  k ,  p ,  and  a  for  the  compact  TEC  model  in  a  prompt  and 
easy  manner.  In  this  concern,  the  convenient  parameter  extraction 
methodology  utilizing  manufacture’s  published  data  [24]  has  been 
widely  extended  in  one  dimensional  (ID)  TEC  modeling  and  design 
[25-27].  The  method  does  not  require  additional  information  ex¬ 
cept  the  readily  available  datasheet  to  produce  the  TEC  perfor¬ 
mance  against  applied  conditions.  However,  this  method  only 
extracts  ID  lumped  parameters  of  the  TEC  module  such  as  thermal 
and  electrical  resistance,  but  does  not  give  the  effective  material 


properties  for  the  3D  compact  model.  So  far,  a  comprehensive 
co-design  methodology  of  compact  modeling  and  parameter 
extraction  of  TECs  in  the  CFD  approach  has  not  been  available. 

This  paper  presents  a  compact  TEC  model  running  on  the  com¬ 
mercially  available  CFD  code  FLUENT,  on  which  the  innumerous 
thermal  management  models  already  existing  can  be  efficiently 
integrated  with  TECs.  The  compact  model  is  accompanied  with 
an  analytical  scheme  of  parameter  extraction,  where  only  four 
parameters  from  TEC  manufacturers’  datasheets  (ATmax,  /max, 
Q.max>  Tit)  are  used  to  characterize  k ,  p ,  and  a  of  the  “black  box”, 
whilst  other  information  of  the  device  is  not  required.  In  order  to 
generalize  the  methodology,  the  detailed  model  of  thermoelectric 
generators  [15]  is  extended  to  be  a  TEC  model  with  full  thermocou¬ 
ple  geometry.  The  physical  TEC  model  is  then  used  to  virtually  ex¬ 
tract  the  effective  k ,  p,  and  a  of  the  compact  model  for  the  general 
scenario  where  all  TEC  parameters  are  completely  available.  The 
physical  model  and  the  compact  model  give  almost  exactly  the 
same  result,  but  in  combination  with  an  air-cooled  heat  exchanger 
as  an  example  analyzing  cooling  enhancement,  the  latter  improves 
the  computation  efficiency  significantly. 


2.  Compact  TEC  model  and  analytical  parameter  extraction 

A  typical  TEC  module  is  shown  in  Fig.  1(a),  in  which  many  n-  and 
p- type  semiconductor  legs  composing  the  cooler  are  connected 
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(b) 


Fig.  1.  TEC  module,  (a)  Detailed  full  geometry,  (b)  Compact  modeling  scheme. 


thermally  in  parallel  between  the  hot  and  cold  ceramic  substrates, 
and  electrically  in  series  by  the  interconnect  metals  to  receive  the 
current  from  the  external  power  supply.  Due  to  the  thermoelectric 
effects,  heat  is  removed  from  the  cold  side  of  the  TEC,  which  is  at¬ 
tached  to  the  object  to  be  cooled,  while  the  TEC  hot  side  dissipates 
heat  to  the  ambient.  The  net  effect  is  that  a  temperature  gradient 
is  created  within  the  TEC. 

Fig.  1(b)  presents  the  compact  scheme  of  the  full  model  in 
Fig.  1(a).  There  are  only  three  sections,  in  which  the  bottom  and 
top  blocks  correspond  to  two  ceramic  substrates,  and  the  macro 
block  in  the  middle  represents  all  thermoelements  and  the  inter¬ 
connect  metal  bridges  above  and  below  the  n-  and  p- type  elements 
of  thermoelectric  material.  Compared  with  using  the  detailed 
geometry,  the  compact  model  is  obviously  advantageous  for  carry¬ 
ing  out  computationally  intensive  parametric  studies  across  vari¬ 
ous  length  scales  during  system  prototype  design.  However,  as 
the  compact  model  greatly  simplifies  the  complicated  TEC  device 
structure  into  a  regular  cuboid  careful  extraction  and  implementa¬ 
tion  of  the  effective  k ,  p,  and  a  are  critical  to  ensure  the  block  to 
faithfully  characterize  the  TEC. 

The  characteristics  of  a  particular  TEC  that  define  its  perfor¬ 
mance  are  ATmax,  Jmax,  and  Qmax  (defined  explicitly  below).  While 
the  thermoelectric  materials  ZT  primarily  determines  A Tmax,  the 
extensive  values  /max,  and  particularly  Qmax  depend  on  the  device 
size,  primarily  through  the  area  A ,  and  height  /,  of  the  TEC.  For  opti¬ 
mum  sizing  of  a  TEC,  it  is  convenient  to  formulate  the  heat  flow 
equations  so  that  the  adjustable  geometric  variables  are  separated 
from  the  internal,  fixed  material  and  device  properties.  In  the  fol¬ 
lowing  analysis,  the  effects  of  the  hot  and  cold  substrates  are  not 
immediately  included  to  directly  utilize  ATmax,  Jmax,  and  Qmax  etc, 
specified  by  the  TEC  manufacturer.  The  effect  of  substrates  will 


be  discussed  at  the  end  of  the  analysis,  as  an  additional  thermal 
conduction  element. 

We  start  with  the  traditional  formulation  for  TEC  by  assuming 
temperature  independent  properties  of  thermoelectric  materials. 
The  rate  of  heat  being  removed  from  the  object  being  cooled  (with 
temperature  Tc)  of  the  TEC  is  written  as, 

Q  =  SITc-KAT-^I2R,  (1) 


where  the  three  terms  at  the  right  side  represent  the  contribution  of 
the  Peltier  effect,  the  heat  conduction,  and  the  Joule  effect,  respec¬ 
tively  [28].  Although  the  temperature  dependence  of  material  prop¬ 
erties  can  be  significant  in  the  evaluation  of  thermoelectric  device 
performance  [15],  the  common  assumption  of  constant  material 
attributes  made  in  compact  TEC  models  usually  does  not  give  rise 
to  a  main  source  of  error  because  the  temperature  gradient  in  TEC 
is  not  so  large  when  compared  with  the  case  of  thermoelectric  gen¬ 
erators  unless  operating  near  ATmax.  When  the  thermal  load  Q  is 
large,  the  side  being  cooled  may  actually  be  at  a  higher  temperature 
than  the  side  being  heated,  making  AT  <  0.  This  is  also  the  case 
when  the  TEC  is  working  as  a  power  generator. 

The  electrical  current  in  the  device  that  gives  the  maximum  AT 
can  be  obtained  by  the  derivative  of  AT  with  respect  to  I  using  (1 ), 


dAT 

~dT 


=  0  =^>  /m 


5Tc 

~T' 


When  Q  =  0, 


(2) 


S2T2 

AT  =  c 
max  2  KR  ' 

Substituting  Jmax  into  (1)  yields 


(3) 


S2T2 

Q  =  ^-KAT.  (4) 

Qmax  of  the  TEC  is  the  Q.  achieved  when  AT  =  0.  It  is  found  to  be 


s2t2 

~2R~' 


(5) 


At  this  point  the  parameters  needed  for  calculation  of  the  per¬ 
formance  of  the  cooler  (S,  I<,  and  R  in  Eq.  (1))  can  be  determined 
from  the  appropriately  defined  specifications  of  cooler  perfor- 


I<  =  and  R  =  2% 

max  l„ 


and  S  =  2  ^ 


Jmax  Th-ATmax 

ATmax  can  be  further  solved  giving  AT, 


in  terms  of  Th  which  is 


more  commonly  controlled  than  Tc  in  a  TEC, 


ATmax  = 


ZTi 


1  +ZTh  +  VI  +2ZTh 


(6) 


Equation  (1)  is  succinctly  rewritten  in  terms  of  the  extensive  vari¬ 
ables  relative  to  their  maximum  quantities  to  obtain  (again  for  con¬ 
stant  Tc  as  opposed  to  the  more  commonly  used  constant  Th) 


J— Lf2— 

Q_max  ^max  \  fmax  /  ATmax 


(7) 


2.1.  Internal  model  of  TEC 

When  the  intensive  TEC  material  properties  are  known  or  de¬ 
sired,  they  can  be  related  to  the  extensive  properties  of  the  TEC 
(S,  I<,  and  R).  S  =  2Not  where  2ot  =  otp-  otn.  Note  that  for  similar  n 
and  p- type  thermoelectric  materials  otn  «  -a p,  and  a  can  then 
approximately  be  the  average  thermopower  of  the  thermoelectric 
materials  (a  =  -  \(xn  «  ap). 
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The  TEC  thermal  conductance,  I<,  can  be  considered  the  sum  of 
all  the  parallel  thermal  conductance  paths  between  the  hot  and 
cold  sides, 

K  =  Nf^nAn  +  KI^L  +  K\  (8) 

Here  is  the  compilation  of  parallel  thermal  losses  associated 
with  the  underfill,  egg-crate,  radiation,  gas  conduction  etc.,  per 
couple.  The  relative  area  ypn  is  fixed  by  optimizing  TEC  perfor¬ 
mance.  The  total  TEC  area,  A ,  depends  not  only  on  the  number  of 
couples,  N  and  the  area  of  the  thermoelectric  elements  An  +  Ap, 
but  also  /,  the  fraction  of  A  covered  by  the  TE  elements  of  the  TEC, 


Af  =  N(A„+AP). 

(9) 

Thus  1<  can  be  simplified  to 

II 

(10) 

where 

k  =  Knyn  +  Kpyp  +  ku 

(11) 

with  Ki  =  %jl<i  again  due  to  parallel  thermal  conduction  losses. 
yp  =  and  yn  =  are  the  relative  areas  for  the  p- type  and 
the  n-type  thermoelectric  elements,  respectively.  Note  that  for  sim¬ 
ilar  n  and  p- type  thermoelectric  materials,  k  is  approximately  the 
thermal  conductivity  of  the  thermoelectric  materials 

(k  =  1/C„-1kp  Ri  Kp). 

Similarly,  ft  is  the  sum  of  the  resistances  of  the  n-  and  p- type 
elements,  and  contact  plus  strap  (interconnect),  ft/,  per  couple, 


Af  a2T2 

T~2p”' 


(19) 


Note  that  unlike  Qmax,  ATmax  does  not  scale  with  any  device 
geometric  parameter  (, A ,  If  or  N),  it  is  only  a  function  of  the  mate¬ 
rials  properties  through  the  thermoelectric  figure  of  merit.  The 
heat  flux  ((ll A)  however  scales  with  the  packing  fraction  and 
length  C/70  and  not  with  the  number  of  elements  n.  The  number 
of  elements  only  becomes  important  when  the  absolute  current 
(as  opposed  to  reduced  current  i)  or  resistance  is  needed.  The  cur¬ 
rent  scales  as  Af/Nl ,  (Eq.  (15))  while  the  resistance  scales  as  N2l/Af. 
The  N 2  scaling  of  resistance  is  necessary  because  a  doubling  of  the 
number  of  thermoelectric  elements  (keeping  area  and  filling  frac¬ 
tion  constant)  requires  also  a  halving  of  the  area  of  each  element, 
which  doubles  the  resistance  per  element. 

The  effective  average  thermoelectric  properties  a,  p,  k ,  can  now 
be  derived  from  manufacturer  A Tmax,  Jmax,  Qmax  specifications  and 
TEC  geometry  ( A  and  /)  using  a  reformulation  of  the  mathematic 
expressions  (15),  (16),  and  (19)  as  long  as  the  filling  fraction /(typ¬ 
ically  about  0.25)  and  number  of  couples  N  are  assumed,  namely: 


a  = 


Qmax 

NTcImax 


(20) 


P  = 


Af  Qm, 
21  N2ll 


(21) 


l  Qmax 

AfATm^ 


(22) 


R  =  n(^  +  ^'  +  R().  (12) 

The  resistance  simplifies  to 
N2l 

R  =  Jf4  P>  (13) 

where 

4  p  =  ^  +  ^  +  R,,  (14) 

in  ip 

with  pl  =  4£ft/.  Note  that  for  similar  n  and  p- type  thermoelectric 
materials,  p  is  approximately  the  effective  resistivity  of  the  ther¬ 
moelectric  materials  (p&\  (j^  +  yi)  ~  Pp)- 
The  /max  of  Eq.  (2)  is 

j  _Af_uh  M5) 

max  Nl  2p  ’  1  j 

and  A Tmax  of  Eq.  (3)  is 

ZT2 

ATmax  =  ^,  (16) 

2  a2 

where  Z  =  %-  approaches  zp  n  =  .  pf  when  there  are  no  electrical  or 

PK  ri  Ppp^P,™ 

thermal  losses,  and  similar  n  and  p- type  thermoelectric  materials 
assuming  temperature  independent  properties. 

Using  the  average  material  parameters,  Eq.  (1)  is  succinctly 
rewritten  in  terms  of  the  reduced  current  i, 


l 

1  ~  i 

lmax 

to  obtain  (for  constant  Tc) 

When  i  =  1,  AT  =  0,  Qmax  is  obtained, 


(17) 


(18) 


However,  the  mathematic  manipulation  above  must  be  cali¬ 
brated  because  of  the  different  measurements  of  Tc  in  ATmax,  Jmax, 
and  Qmax-  As  is  stated,  the  condition  in  determining  Qmax  in  (5)  and 
(19)  is  AT  =  0,  i.e.,  Th  =  Tc,  whilst  the  objective  of  Jmax  is  to  maximize 
AT,  i.e.,  to  achieve  ATmax.  Since  TEC  performance  data  are  com¬ 
monly  published  for  constant  hot  side  temperatures  (usually 
27  °C  and  50  °C),  (19),  e.g.,  is  changed  to 


Af  oc2Tl 

l  2  p 


(23) 


in  order  to  conform  with  the  expressions  of  Jmax  and  Qmax, 
where  Tc  =  Th-  A Tmax.  Combining  (15),  (16),  and  (23)  yields 


a  - 


Qmax(Th  ATit 


(24) 


P  = 


Af(Jh  -  AT, 


2  T2hl 


N  It 


(25) 


7/.  _  K^h  ~  ATmax)  Q.max 

K~~ml  ^  11 

In  (24)-(26),  ATmax,  /max,  (lmax,  and  Th  are  directly  available  from 
the  TEC  vendor.  The  number  of  elements  and  element  dimensions 
can  also  be  easily  derived  from  the  device  diagram  in  datasheets  or 
measured  from  the  real  module.  Thus,  the  effective  material  prop¬ 
erties  a,  p,  k,  for  the  compact  model  are  successfully  extracted. 

In  the  above  formulation,  parallel  heat  loses  due  to  I<i  can  sim¬ 
ply  be  included  by  increasing  k.  In  a  real  TEC  there  is  also  thermal 
resistance  between  the  external  surfaces  of  the  TEC  and  the  inter¬ 
nal  junction  between  the  thermoelectric  materials  and  the  inter¬ 
connect  metals.  As  is  shown  in  Fig.  2,  all  the  series  thermal 
resistances  between  the  external  cold  temperature  measurement 
(Tc)  and  the  internal  thermoelectric  material  cold  junction  temper¬ 
ature  (Td)  are  summed  together  to  give  <9C.  This  may  include  ther¬ 
mal  resistances  from  the  substrate,  interconnects,  thermal  contact 
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Fig.  2.  1 D  heat  transfer  analogy  from  the  heat  source  to  the  heat  sink. 


resistance  and  substrate  spreading  resistance.  Similarly  defined  on 
the  hot  side  are  Th,  Thi,  and  © h . 

The  importance  of  these  thermal  resistances  depends  on  their 
relative  magnitude  compared  to  the  thermal  resistance  of  the  ther¬ 
moelectric  material.  Thus  these  additional  thermal  resistances  be¬ 
come  more  important  as  the  height  of  the  thermoelectric  elements 
is  decreased.  However,  ©c  and  0 h  are  usually  unknown  in  manu¬ 
facturers’  datasheets,  so  the  maximum  of  the  internal  temperature 
difference,  A Tte,  cannot  be  obtained  from  the  published  external 
A Tmax.  Moreover,  the  current  value  Jmax  which  maximizes  Th  -  Tc, 
does  not  necessarily  maximize  Thi  -  Tci. 

We  notice  that  A Tmax,  Jmax,  Q.max  and  Th  in  datasheets  are  exter¬ 
nal  performance  parameters  measured  for  the  entire  TEC,  where 
the  effects  of  ©c  and  ©h  are  included.  When  one  is  evaluating 
TEC  through  a  compact  model,  the  main  concern  is  to  evaluate 
parameters  that  reproduce  device  performance  accurately.  There¬ 
fore,  in  this  compact  model,  the  effective  a,  p,  k  of  the  main  block 
containing  the  thermoelectric  materials  as  opposed  to  the  top/bot¬ 
tom  layers  are  directly  calculated  by  manufacturers’  A Tmax,  Jmax, 
Qmax  and  Th  using  (24)-(26).  In  this  way,  all  thermal  and  electrical 
properties  of  the  TEC,  reflected  by  the  performance  parameters,  are 
assumed  to  result  entirely  from  the  thermoelectric  layer.  To  test 
these  values  in  the  model,  a  very  large  thermal  conductivity,  much 
larger  than  the  ceramic  material,  is  applied  to  the  top  and  bottom 
substrates  of  the  test  structure  to  essentially  eliminate  the  effect  of 
the  thermal  resistances  © c  and  © h .  Because  very  large  values  of 
conductivity  inhibit  rapid  convergence  of  the  calculation,  a  value 
of  200  W  m_1 1<-1  is  applied  to  the  two  substrates  in  the  test  cases. 

The  3D  TEC  model  is  implemented  in  the  64  bit  ANSYS  FLUENT 
12.1.  User-defined  functions  (UDFs)  are  used  to  model  the  effective 
heat  conduction,  volumetric  Joule  heat  source,  and  interfacial  Pel- 
tier  heat  in  terms  of  a,  p,  k  extracted  from  A Tmax,  Jmax,  Qmax  and  Th 
in  datasheets.  The  UDF  for  the  effective  conductivity  of  the  TEC 
materials  is  implemented  in  terms  of  (26)  to  model  the  Fourier 
heat  flow  in  the  middle  TE  block. 

An  energy_source  UDF  is  written  to  uniformly  generate  the  vol¬ 
umetric  Joule  heat  on  all  cells  of  the  middle  block.  By  combing  (13) 
and  (25),  the  Joule  heating  density  can  be  obtained  as 

_I2R_  I2  N24pl  _  2  2 (Th  -  A Tmax)2Qmax 
q>oul  Al  Al  Af  mj\  '  K  ] 

where  7(i)  can  be  determined  as  the  model  input  parameter  repre¬ 
senting  the  electrical  current  applied  to  the  TEC. 

The  Peltier  heat  is  absorbed  and  evolved  at  the  cold  and  hot 
junctions  of  dissimilar  materials,  respectively.  In  the  compact 
model,  two  Peltier  UDFs  are  used  to  represent  the  absorbed  and 
evolved  heat,  and  placed  at  the  two  interfaces  between  the  middle 


1(A) 


Fig.  3.  TEC  temperature  difference  comparison  between  compact  simulation  and 
datasheet  for  various  currents  (a)  Q.  =  0  W,  Th  =  300  K;  (b)  Q.  =  20  W,  Th  =  300  K. 


0.1  0.7  1.3  1.9  2.5  3.1  3.7  4.3 

1(A) 


Fig.  4.  TEC  power  supply  voltage  comparison  between  compact  simulation  and 
datasheet  for  various  currents  (Q.  =  0  W,  Th  =  300  K). 

TE  block  and  the  bottom/top  substrates  as  their  boundary  profiles. 
FLUENT  is  able  to  solve  a  ID  conduction  equation  of  the  interfacial 
walls  as  well  as  the  heat  generation  in  the  wall.  To  include  Peltier 
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Fig.  5.  Temperature  contours  of  (a)  3D  physical  TEC  model,  (b)  3D  compact  TEC  model,  (c)  2D  cold  side  by  physical  TEC  model.  (Unit:  K,  Q.  =  20  W,  Th  =  300  K,  /  =  5.7  A). 


Table  1 

Summary  of  the  physical  TEC  model. 


Zone  type 

Material 

Source  terms 

Potential 

fields 

Solid 

number 

p- Type  leg 

p-Type  thermoelectric  Joule;  Thomson;  Seebeck; 

127 

material 

Peltier 

Ohm 

n-Type  leg 

n-Type  thermoelectric  Joule;  Thomson;  Seebeck; 

127 

material 

Peltier 

Ohm 

Top  bridge 

Conducting  metal 

Joule 

Ohm 

127 

Bottom 

Conducting  metal 

Joule 

Ohm 

126 

bridge 

Bottom 

Conducting  metal 

Joule 

Ohm 

1 

bridge  (+) 

Bottom 

Conducting  metal 

Joule 

Ohm 

1 

bridge  (-) 

Substrate 

Ceramic 

2 

effects  in  the  heat  transfer  calculation,  the  thickness  of  the  wall  is 
specified  as  10-5  m.  We  note  that  a  non-zero  thickness  must  be 
specified  in  FLUENT  for  such  thin  wall  thermal  calculation,  even 
though  this  thickness  does  not  really  exist  in  the  solid  model.  By 


using  (24),  the  heat  generation  rate  due  to  Peltier  effect  can  be  ob¬ 
tained  as 

N2CClThi,ci  Qmax(Th  -  AT  max  )  N2/rw,d 

%elt  ±  a,1()_5  ±  m2Imax  a.1q-5 

_  Q.max  (T„  -  AT  max  )27yd 

t\k io  5 

where  the  plus  and  minus  signs  stipulate  whether  the  Peltier  heat  is 
evolved  or  absorbed.  The  UDFs  are  able  to  identify  the  interfaces 
and  external  surfaces  automatically,  and  hence  calculate  the  aver¬ 
age  Th,  Tc,  Thi  and  Tci  at  the  four  surfaces  as  the  simulation  results. 
Both  the  internal  A Tte  and  the  external  AT  can  be  similarly  calcu¬ 
lated  at  the  end  of  each  iteration,  and  V  is  also  obtained  in  terms 
of  the  internal  A Tte, 

V  =  SATte  +  RI.  (29) 

In  order  to  validate  the  compact  model  as  well  as  the  parameter 
extraction  method,  the  TEC  DTI  2-4  made  by  Marlow  is  simulated 
as  the  example.  The  TEC  has  127  thermocouples  (N=  127)  and  a 
30  *  30  mm2  cross-section,  and  hence  a  filling  factor  of  0.2822. 
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The  four  key  performance  parameters  are  described  in  the  manu¬ 
facturer’s  datasheet  when  the  hot  temperature  is  fixed  at  300  K, 
i.e.,  ATmax  =  66  K,  /max  =  3.7  A,  Q.max  =  36  W  [29].  Substituting  these 
values  into  (24)-(26),  effective  material  properties  for  the  compact 
model  can  be  obtained,  i.e.,  a  =  199pV/K,  p  =  13pQm, 
k=  1.31  Witt1  1C1. 

Fig.  3  and  4  compare  performance  curves  published  in  the  data¬ 
sheet  [29]  and  the  simulation  results  of  the  compact  model  for 
Th  =  300  K.  The  scenario  of  zero  thermal  load  (Q.  =  0  W)  is  presented 
in  Fig.  3(a),  where  AT  changes  with  J.  Fig.  3(b)  presents  the  case 
that  Q  =  20  W,  which  is  closer  to  the  realistic  cooling  applications. 
In  Fig.  4,  the  change  of  V  with  /  is  depicted  for  comparison  between 
the  datasheet  and  the  simulation.  For  most  current  values,  it  can  be 
seen  that  the  simulated  AT  and  V  match  the  datasheet’s  character¬ 
istics  very  well.  The  deviation  becomes  slightly  obvious  when  Q  is 
increased  and  when  I  passes  its  maximum  value  (J  >  Jmax),  likely 
caused  by  measurement  uncertainties  of  performance  parameters 
and  the  simplification  of  the  compact  modeling  approach.  How¬ 
ever,  given  the  benefits  on  computational  cost  brought  by  the  com¬ 
pact  model,  such  a  small  error  level  is  fairly  acceptable  for  practical 
cooling  system  design  [16-23]. 


instead  of  (29),  where  the  first  term  is  the  sum  of  Seebeck  po¬ 
tential  of  all  thermoelements  working  against  /,  and  the  second 
term  is  the  ohmic  voltage  drop  as  /  flows  through  all  thermoele¬ 
ments  of  the  TEC  module.  aPi„  and  pp>n  are  dependent  properties 
of  the  3D  temperature  distribution  T  for  each  n-  or  p- type  leg  in 
the  physical  model.  The  starting  and  the  end  points  of  the  curvilin¬ 
ear  integral  F  depend  on  the  actual  contacts  in  the  junctions  be¬ 
tween  which  the  effective  Seebeck  potential  and  resistance  of 
each  leg  are  measured.  In  this  model  the  voltage  drop  across  the 
strap  (Vstrap)  resistance  of  the  bridges  are  also  included  by  the  con¬ 
tinuous  potential  field  although  not  described  in  detail  in  (30). 

Similarly,  the  full  3D  temperature  field  of  the  TEC  can  be  estab¬ 
lished  by  the  physical  model,  as  is  shown  in  Fig.  5(a),  in  which  the 
/-related  source  terms  of  Peltier,  Thomson,  and  Joule  heat  are  all  in¬ 
cluded.  The  TEC  thermal  load  Q  in  the  physical  model  consists  of 
the  absorbed  Peltier  heat  subtracted  by  the  nonlinearly  conducted 
heat  from  the  hot  side,  i.e., 


Q  =  E  /// Vapr;dn+ J j J I -KpVTdAp- J j -K„VTdA„ 


(31) 


3.  Physical  TEC  model  and  virtual  parameter  extraction 

To  determine  quantitatively  how  much  computational  cost  can 
be  saved  by  using  the  compact  model,  it  is  necessary  to  compare 
the  same  TEC  in  a  modeling  approach  that  uses  the  full  geometry. 
The  3D  thermoelectric  model  in  [15],  validated  by  experimental  re¬ 
sults,  is  ideal  for  this  purpose  because  it  is  also  implemented  in 
FLUENT,  so  that  any  comparisons  should  be  computationally  fair. 
Also,  there  may  be  a  need  to  extract  the  effective  material  proper¬ 
ties  from  the  measured  device  properties  used  in  the  compact 
model,  especially  for  those  TECs  under  laboratory  development 
and  testing.  Although  manufacturer’s  specifications  are  a  reliable 
data  source,  a  full  TEC  model  is  important  to  evaluate  also  in  the 
case  when  all  material  parameters,  including  the  temperature 
dependency  of  a,  p,  k ,  thermal  and  electrical  contact  resistance 
are  given  instead  of  the  device  performance.  The  full  TEC  model 
can  be  used  to  calculate  ATmax,  Jmax,  Q.max  for  a  given  Th  from  the 
material  properties  and  geometry. 

As  the  3D  model  in  [15]  faithfully  includes  all  coupled  thermo¬ 
electric  as  well  as  components  that  provide  losses  and  other  para¬ 
sitic  effects,  it  can  be  referred  to  as  the  physical  model  in  contrast 
with  the  compact  model  which  is  only  concerned  with  overall  per¬ 
formance.  The  physical  model  used  here  was  originally  created  for 
thermoelectric  generators  (TEG)  [15],  so  the  model  should  be 
transplanted  for  coolers  as  the  first  step.  Due  to  the  same  govern¬ 
ing  equations  of  thermoelectricity  for  TEG  and  TEC,  the  numerical 
schemes  and  model  framework  are  the  same  except  that  the  con¬ 
cept  of  a  load  resistance  of  a  TEG  is  replaced  by  the  the  current  sup¬ 
plied  I  used  for  TECs. 

The  3D  physical  model  of  a  127  couple  (N=  127)  TEC  module  is 
shown  in  Fig.  5(a)  and  the  modeling  domain  is  summarized  in  Ta¬ 
ble  1.  The  minus  and  plus  pads,  which  used  to  simulate  the  lead 
wires  of  the  TEG,  are  now  used  to  pass  through  the  prescribed  I 
to  establish  the  electric  current  density  through  the  TEC.  The  cur¬ 
rent  direction  is  applied  in  a  way  to  ensure  that  Peltier  heat  is  ab¬ 
sorbed  at  the  cooling  side  and  released  at  the  other. 

In  the  physical  TEC  model,  V  of  the  external  power  supply  is  cal¬ 
culated  as 

V=Y,(yJ  UpVTdr  -  j  oinVTdr'j 

+%(jTpvwxpdr+STp^kdr)+Vstmr'  (30) 


where  ap>n  and  Kp<ri  are  dependent  properties  of  T  for  each  n-  or  p- 
type  leg.  The  volume  integral  Q  and  the  surface  integral  Ap<n  are  with 
respect  to  Peltier  heat  density  and  the  conduction  heat  flux  at  the 
n-  or  p- type  cold  junctions. 


1(A) 


Fig.  6.  TEC  temperature  difference  comparison  between  compact  simulation  and 
full  simulation  for  various  currents  (a)  Q.  =  0  W,  Th  =  300  K.(b)  Q.  =  20  W,  Th  =  300 IC. 
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Within  FLUENT,  the  3D  physical  model  solves  the  nonlinearly 
coupled  thermal  and  electrical  processes  such  as  (30)-(31)  for  the 
components  in  Tabel  1.  The  experimental  results  used  to  deter¬ 
mine  all  the  TEC  performance  parameters  can  be  measured  virtu¬ 
ally.  The  thermoelectric  module  TEC1 -12706  made  by  Tianjin 
Institute  of  Power  Sources,  China,  is  chosen  in  this  section  as  the 
simulation  example  because  the  detailed  temperature  dependenc¬ 
es  of  the  n-  and  p- type  bismuth  telluride  materials  are  available 
from  the  manufacturer.  The  thermal  and  electrical  contact  resis¬ 
tances  of  TEC1  -12706  are  also  experimentally  determined  in  the 
previous  study  [30]. 

Unless  a  uniform  distribution  is  imposed,  the  temperature  pro¬ 
file  of  TEC  ceramic  surfaces  is  usually  non-uniform  primarily  due  to 
the  space  between  TE  elements  within  the  module,  as  is  shown  in 
Fig.  5(c).  Therefore,  in  dealing  with  the  AT  information  associated 
with  the  physical  model,  the  arithmetic  average  of  the  tempera¬ 
tures  of  all  cell  faces  at  the  TEC  surface,  calculated  automatically 
in  the  UDF,  is  actually  plotted  in  the  following  figures. 


The  squares  in  Fig.  6(a)  present  the  physical  simulation  of  the 
change  of  AT  with  I  when  Q  =  0  W.  It  can  be  seen  that 
ATmax  =  57.63  K  occurs  at  the  current  Jmax  =  5.7  A.  Qmax  =  51.25  W 
can  be  obtained  by  another  simulation  in  which  Jmax  is  applied 
and  the  both  sides  of  the  TEC  prescribed  to  be  the  same  tempera¬ 
ture.  In  terms  of  (24)-(26),  the  effective  a,  p}  k  can  be  obtained  for 
the  compact  TEC  model,  as  is  shown  in  Fig.  5(b). 

Fig.  6(a)  and  (b)  show  clearly  that  the  physical  model  and  the 
compact  model  depict  the  AT  change  as  a  function  of  I  with  almost 
the  same  characteristics.  The  small  difference  between  the  two 
simulations  is  slightly  enlarged  when  the  thermal  load  is  increased, 
similar  to  the  results  presented  in  Section  2.  In  Fig.  7(a)  and  (b),  the 
changes  of  V  with  I  by  the  physical  model  and  the  compact  model 
are  compared  under  different  thermal  loads.  For  most  current  val¬ 
ues,  it  can  be  seen  that  the  two  models  match  very  well.  When  Q  is 
increased  and  when  I  passes  its  maximum  value,  once  again,  the 
deviation  becomes  more  obvious  as  the  results  in  Section  2.  It  is 
noteworthy  that,  when  /  >  Jmax,  the  V  -  I  curves  from  the  physical 
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Fig.  7.  TEC  power  supply  voltage  comparison  between  compact  simulation  and  full  simulation  for  various  currents  (a)  Q.  =  0  W,  Th  =  300  K;  (b)  Q  =  20  W,  Th  =  300  K. 
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Fig.  8.  TEC  based  heat  sink  of  a  cooling  system  using  the  compact  CFD  model,  (a)  schematic  structure  and  dimension,  (b)  3D  temperature  distribution,  (c)  2D  temperature 
distribution  of  the  TEC  hot  side  surface,  (Q.  =  0  W,  I  =  5.7  A,  velocity:  30  m/s). 


model  shows  a  nonlinear  change,  which  cannot  be  accurately  cap¬ 
tured  by  the  compact  model.  As  explained  above,  the  physical 
model  faithfully  reflects  the  temperature  dependent  material  prop¬ 
erties  and  the  full  coupled  multi-physics  effects.  The  relative  devi¬ 
ations  between  the  two  models  are  also  given  in  Fig.  7,  and  it  can 
be  seen  that  the  relative  deviations  are  actually  within  the  normal 
range  when  the  absolute  deviation  is  enlarged  at  extremely  high 
currents. 

The  temperature  variation  within  the  entire  surface  by  the  com¬ 
pact  model  is  very  small  (265.91  ±  0.01  K),  as  shown  in  Fig.  5(b).  In 
Fig.  5(c)  where  the  local  range  of  the  2D  TEC  cold  side  is  applied  to 
the  temperature  scale,  the  physical  model  clearly  displays  its 
excellent  ability  to  accurately  capture  the  horizontal  temperature 
gradient  caused  by  the  heat  spread  [31],  which  is  very  useful  in 
analyzing  spot  cooling  applications.  Therefore,  in  situations  when 
the  computation  accuracy  of  unusual  situations  is  highly  de¬ 
manded  the  physical  model  is  likely  to  be  more  accurate  than 
the  compact  model. 

On  the  other  hand,  the  compact  model  can  display  its  strength 
when  computational  cost  and  efficiency  are  taken  into  account. 
Regarding  modeling  size,  the  FLUENT  computational  domain  of 
the  physical  TEC  model  contains  658,495  volumes  and  881,900 
nodes,  whilst  the  compact  TEC  model  has  121,600  volumes  and 
143,702  nodes.  In  order  to  speed  up  the  convergence  of  the  poten¬ 
tial  field,  double-precision  solver  (3DDP)  is  used  for  the  physical 
model  instead  of  a  single  precision  solver,  and  on  a  Lenovo  W520 
laptop  workstation,  the  model  takes  approximately  1 0  h  in  average 
to  finish  the  computation.  However,  the  compact  TEC  model  takes 
only  a  few  minutes  to  converge  on  the  same  computer  when  3DDP 
is  selected  due  to  the  simplification  on  both  gridding  and  comput¬ 
ing.  Thus,  in  many  scenarios  of  cooling  system  optimization  where 


the  computation  efficiency  is  an  important  factor,  a  sensible  trade¬ 
off  between  the  accuracy  of  the  physical  model  and  the  speed  of 
the  compact  model  should  lead  to  selection  of  the  compact  model 
given  their  small  simulation  differences  under  normal  applications 
as  shown  in  Fig.  6  and  7. 

4.  Model  exploration  at  the  system-level 

While  the  processes  and  phenomena  within  a  TEC  are  complex, 
its  compact  model  has  a  straightforward  and  feasible  path  to  ther¬ 
mal  management  applications  when  all  the  important  effective 
parameters  are  appropriately  extracted.  To  provide  a  clear  descrip¬ 
tion  how  the  compact  model  is  applied  to  the  design  of  practical 
cooling  systems,  an  assembly  comprising  a  TEC  and  an  air  cooled 
heat-sink  is  modeled  in  FLUENT.  As  is  shown  in  Fig.  8(a),  the  TEC 
is  assumed  mounted  on  a  component,  and  the  bottom  surface  of 
the  TEC  represents  the  heat  source  on  which  either  a  heat  flux  or 
a  temperature  can  be  specified.  It  is  assumed  that  air  flows  through 
the  duct  and  fins  of  the  heat  sink  so  that  it  cools  the  TEC  by  forced 
air  convection. 

In  the  CFD  model,  the  finned  geometry  of  the  heat  sink,  as  an 
empirical  design,  indicates  the  presence  of  laminar  airflow.  A  con¬ 
stant  inlet  temperature  of  300  K  and  a  number  of  air  velocities 
were  considered  as  fluid  dynamic  boundary  conditions.  Thermal 
boundary  conditions  on  other  walls  that  are  not  attached  to  the 
heat  source  or  exposed  to  the  airflow  are  assumed  adiabatic.  In 
the  simulation  of  such  a  cooling  system,  the  inputs  are  the  inlet 
air  velocity  and  temperature,  the  source  temperature/heat  genera¬ 
tion  rate,  the  electrical  current  to  the  TEC,  and  the  outputs  are  the 
3D  temperature  and  flow  fields  of  the  whole  computation  domain, 
as  is  shown  in  Fig.  8(b). 
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Fig.  9.  TEC  cold  side  temperatures  (a)  and  thermal  loading  ability  (b)  for  various 
airflow  velocities  (J  =  5.7  A). 

In  the  3D  temperature  contour,  it  can  be  seen  that  there  is  an 
obvious  temperature  gradient  from  the  flow  (outlet)  to  the  heat 
sink  base  due  to  the  small  thermal  resistance  of  the  heat  sink. 
The  TEC  enables  an  “active”  heat  sink  by  reducing  the  cooling  tem¬ 
perature  below  the  ambient  (300  K).  The  effect  of  the  air  flow  on 
the  2D  temperature  distribution  of  the  base  surface  in  the  heat  sink 
is  shown  in  Fig.  8(c),  which  displays  a  non-uniform  profile. 

Fig.  9(a)  shows  how  the  temperature  of  the  TEC  cold  side  can  be 
reduced  by  increasing  the  airflow  velocity  to  enhance  the  convec¬ 
tion.  It  can  be  seen  that  the  optimal  velocity  is  around  10  m/s  for 
various  thermal  loads,  and  when  the  velocity  is  higher  than 
15  m/s,  further  enhancement  of  the  convection  cannot  significantly 
improve  the  cooling  performance.  Similarly,  the  heat  pumping 
ability  from  the  cooling  objective  to  the  flow  fluid  of  the  heat  sink 
is  optimized  within  the  range  5-15  m/s,  after  which  the  thermal 
load  cannot  be  significantly  increased,  as  is  shown  in  Fig.  9(b). 

5.  Conclusions 

This  paper  presents  a  compact  3D  TEC  model  by  which  the  dif¬ 
ferent  materials  and  components  within  the  TEC  module  are  ab¬ 
stracted  as  a  single  block.  A  significant  amount  of  grid  is  reduced 
in  modeling  a  TEC,  and  also  the  grid  size  mismatch  between  the 
TEC  and  thermal  system  is  mitigated.  It  is  also  shown  how  the 
effective  parameters  of  the  compact  model  can  be  extracted  from 
manufacturer’s  data  or  calculated  from  materials  parameters.  The 
compact  model  and  parameter  extraction  are  validated  by  compar¬ 
ing  the  simulation  results  and  the  manufacturer’s  performance 
data  at  various  electrical  and  thermal  conditions  with  excellent 
agreement. 


A  full  TEC  model  is  also  implemented  in  FLUENT,  where  the  ex¬ 
act  device  geometry,  coupled  thermoelectric  effects,  and  the  ther¬ 
mal  gradient  and  electric  potential  are  fully  represented  in  3 
dimensions.  The  physical  model  can  be  used  to  carry  out  virtual 
experiments  to  obtain  TEC  performance  data  to  extract  the  effec¬ 
tive  parameters  for  the  compact  model.  Simulation  comparisons 
validate  that  the  compact  modeling  and  parameter  extraction  yield 
results  almost  as  accurate  as  the  physical  model,  but  with  compu¬ 
tational  speed  roughly  100  times  faster. 

The  compact  model  uses  the  CFD  approach,  so  that  the  integra¬ 
tion  of  TEC  into  many  thermal  management  applications  can  be 
readily  implemented.  A  demonstrative  example  of  successful  utili¬ 
zation  of  the  compact  model  is  presented  by  modeling  an  air¬ 
cooled  heat  sink.  Rapid  optimization  of  the  cooling  system,  e.g., 
fan  arrangement,  airflow  speed  and  angle,  fin  and  spreader  geom¬ 
etry,  as  well  as  the  TEC  operating  current,  can  be  done  in  a  system¬ 
atic  way.  A  detailed  comparison  of  various  spacing,  length,  and 
width  of  the  fins  to  enhance  the  heat  transfer  for  TEC  applications, 
taking  advantage  of  the  proposed  co-simulation  technique,  can  be 
found  in  [32]. 
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